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X-RAY CT SYSTEM AND BEAM-HARDENING POST-PROCESSING METHOD 

BACKGROUND OF THE INVENTION 
5 The present invention relates to a beam-hardening post-processing method for 

correcting the intensity of X-rays to be transmitted by a subject on the basis of phantom 
data, and to an X-ray computed tomography (CT) system. 

An X-ray source employed in an X-ray CT system generates X-rays that fall within 
a certain energy range. On the other hand, an absorption coefficient for X-rays transmitted 

10 by a subject depends on the energy in the X-rays. The larger a length in a subject over 
which X-rays are transmitted, the larger the average energy. This phenomenon is referred 
to as a beam-hardening effect. Consequently, a proportional relationship is not established 
between the intensity of transmitted X-rays, that is, a projection information value and the 
length over which X-rays are transmitted, but a non-linear relationship is. 

15 The beam-hardening effect causes the cupping effect signifying that the intensity 

detected in the center of a reconstructed image is low. A correction coefficient to be used 
to correct projection information values based on which a reconstructed image is produced 
to exhibit a uniform intensity is calculated in relation to each of the channels of the X-ray 
detector, whereby the correction is achieved (refer to, for example, Patent Document 1). 

20 A plurality of cylindrical phantoms having different diameters that are large enough 

to generally cover the entire field of view (FOV) (scan field) defined in the center of an X- 
ray field is scanned for the purpose of higher-precision correction. Projection information 
acquired from the phantoms is used to improve the precision in correction. 
[Patent Document 1] 

25 Japanese Unexamined Patent Publication No. Hei 5(1 993)- 130987 (p.2 and 3, 

Fig. 1 and Fig. 2) 
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By the way, for highly precise correction of projection information values, the 
largest possible number of different projection information values is needed in relation to 
each of the channels of an X-ray detector. Therefore, lots of phantoms that have different 
diameters must be scanned. 
5 On the other hand, the larger the diameter of a phantom, the lower a signal-to-noise 

ratio reflected in an acquired projection information value. If a correction coefficient 
calculated from projection information values reflecting lowered signal-to-noise ratios is 
used to correct projection information, a reconstruction image is like to suffer degradation 
in image quality deriving from a ring artifact or the like. Consequently, when a subject 
10 having large dimensions is scanned, image quality is likely to be degraded. 

SUMMARY OF THE INVENTION 

An object of the present invention is to provide an X-ray CT system in which even 
when a subject having relatively large dimensions is scanned, projection data can be 

15 highly precisely corrected in terms of the beam-hardening effect in relation to each of the 
channels of an X-ray detector, and to provide a beam-hardening post-processing method 
for X-ray CT systems. 

According to the present invention, there is provided an X-ray CT system in which: 
an X-ray detector having multiple channels is used to acquire projection information, 

20 which is provided by an X-ray beam passing through a scan field, as a plurality of views 
from plural directions; and projection information detected on each of the channels is 
corrected in terms of the beam-hardening effect. The X-ray CT system comprises: a 
correction coefficient producing means that calculates a correction coefficient, which is 
used for correction, from projection information acquired from a phantom placed in the 

25 scan field; a correction coefficient modifying means that uses a first correction coefficient, 
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which the correction coefficient producing means calculates using projection information 
acquired from a first phantom, to modify a second correction coefficient which the 
correction coefficient producing means calculates using projection information acquired 
from a second phantom larger in dimensions than the first phantom; and a correcting 
means that corrects projection information, which is acquired from a subject positioned in 
the scan field, using the first correction coefficient and the corrected second correction 
coefficient. 

Preferably, the correction coefficient producing means comprises: a producing 
means that samples first projection information from projection information acquired from 
a phantom in relation to all views so as to produce one sinogram; a beam-hardening 
correction means that corrects the first projection information in terms of the beam- 
hardening effect so as to produce second projection information; a first fitting means that 
fits a first function to the second projection information so as to produce third projection 
information; and a second fitting means that fits a second function to the third projection 
information values which are provided as functions having as independent variables the 
second projection information values sampled in relation to all the views and each of the 
channels of the X-ray detector. 

More preferably, the X-ray detector is formed with a plurality of detection modules 
each having a predetermined number of channels. The correction coefficient modifying 
means separates the reflections of high-frequency components from correction coefficient 
data calculated using the second phantom so as to leave the dependencies on the detection 
characteristics of the detection modules. The correction coefficient modifying means then 
synthesizes the correction coefficient data, from which the reflections of high-frequency 
components are separated, with the reflections of high-frequency components in correction 
coefficient data calculated using the first phantom. 
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A beam-hardening post-processing method for X-ray CT systems in accordance 
with the present invention comprises: an acquiring step of scanning a first phantom and a 
second phantom, which is larger in dimensions than the first phantom, placed between an 
X-ray tube and an X-ray detector, and acquiring projection information as a plurality of 
5 views from plural directions using an X-ray detector that has multiple channels; a 
producing step of calculating first and second correction coefficients, which are used to 
correct the projection information detected on each of the channels of the X-ray detector, 
from projection information acquired from the first and second phantoms; a modifying 
step of modifying the second correction coefficient using the first correction coefficient; 

10 and a correcting step of correcting projection information, which is acquired from a 
subject positioned in the scan field, using the first correction coefficient and the modified 
second correction coefficient. 

According to the present invention, the first phantom and the second phantom 
larger in dimensions than the first phantom are scanned in order to acquire projection 

15 information. The first and second correction coefficients that are used for correction are 
calculated from the projection information. 

The larger the dimensions of a phantom, the lower a signal-to-noise ratio. 
Therefore, a CT image produced by performing correction using the second correction 
coefficient is more likely to suffer degradation in image quality attributable to a ring 

20 artifacts or the like than a CT image produced by performing correction using the first 
correction coefficient. 

According to the present invention, the first correction coefficient calculated from 
projection information that reflects relatively high signal-to-noise ratios is used to modify 
the second correction coefficient. 

25 Consequently, even when a correction coefficient is calculated using a phantom 
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that has large dimensions, correction can be achieved highly precisely. 

According to the present invention, the precision in correcting projection 
information acquired from a large subject can be improved. 

Further objects and advantages of the present invention will be apparent from the 
following description of the preferred embodiments of the invention as illustrated in the 
accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram showing the overall configuration of an X-ray CT system. 

Fig. 2 is a functional block diagram showing a data processing unit. 

Fig. 3 is a flowchart describing actions to be performed in the data processing unit. 

Fig. 4 shows a sinogram produced from projection information values acquired 
from a phantom, and the projection information values. 

Fig. 5 shows the positional relationship between a phantom and a rotary assembly. 

Fig. 6 is a block diagram showing files stored in a storage device. 

Fig. 7 presents processing to be performed on projection information values in the 
direction of channels. 

Fig. 8 presents processing to be performed on projection information values in the 
direction of views. 

Fig. 9 shows an example of a correction coefficient calculated from projection 
information values. 

Fig. 10 shows examples of correction coefficients calculated from projection 
information values acquired from phantoms having different diameters. 

Fig. 11 includes graphs for explaining an example of a method of modifying 
correction coefficient information. 
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Fig. 12 is a graph presenting an example of correction coefficients calculated using 
two phantoms having different diameters. 

Fig. 1 3 is a graph presenting an example of fitting of a high-order function. 

5 DETAILED DESCRIPTION OF THE INVENTION 

Referring to appended drawings, a preferred embodiment of a beam-hardening 
post-processing method and an X-ray CT system in accordance with the present invention 
will be described below. 

To begin with, the overall configuration of an X-ray CT system of an embodiment 
10 will be described below. Fig. 1 is a block diagram of an X-ray CT system. As shown in 
Fig. 1, the X-ray CT system comprises a scanner gantry 2 and an operating console 6. 

The scanner gantry 2 includes an X-ray tube 20. X-rays radiated from the X-ray 
tube 20 and not shown are recomposed into, for example, a fan-shaped X-ray beam, that 
is, fan-beam X-rays by means of a collimator 22, and then irradiated to an X-ray detector 
15 24. 

The X-ray detector 24 has a plurality of X-ray detection elements set in array in a 
direction in which the fan-beam X-rays spread. The X-ray detector 24 is therefore a multi- 
channel detector having the plurality of X-ray detection elements set in array. 

The X-ray detector 24 forms an X-ray incidence surface curved like a cylindrical 

20 concave surface as a whole. The X-ray detector 24 is formed using, for example, a 
combination of a scintillator and a photodiode. The present invention is not limited to this 
combination. For example, a semiconductor X-ray detection element that utilizes 
cadmium telluride (CdTe) or an ion-chamber type X-ray detection element that utilizes 
xenon gas will do. The X-ray tube 20, collimator 22, and X-ray detector 24 constitute an 

25 X-ray irradiation/detection unit. 
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A data acquisition unit 26 is connected to the X-ray detector 24. The data 
acquisition unit 26 acquires data items detected by the respective X-ray detection elements 
constituting the X-ray detector 24. An X-ray controller 28 controls X-irradiation from the 
X-ray tube 20. Illustrating the connective relationship between the X-ray tube 20 and X- 
5 ray controller 28 and the connective relationship between the collimator 22 and a 
collimator controller 30 will be omitted. The collimator controller 30 controls the 
collimator 22. 

The foregoing components started with the X-ray tube 20 and ended with the 
collimator controller 30 are incorporated in a rotary assembly 34 of the scanner gantry 2. 

10 A subject or a phantom is mounted on a cradle within a bore 29 located in the center of the 
rotary assembly 34. The rotary assembly 34 is rotated while being controlled by a rotation 
controller 36. X-rays are irradiated from the X-ray tube 20. The X-ray detector 24 detects 
X-rays transmitted by the subject or phantom as projection information composed of 
views. Illustrating the connective relationship between the rotary assembly 34 and 

1 5 rotation controller 36 will be omitted. 

The operating console 6 includes a data processing unit 60. The data processing 
unit 60 is realized with, for example, a computer. A control interface 62 is connected to 
the data processing unit 60, and the scanner gantry 2 is connected to the control interface 
62. The data processing unit 60 controls the scanner gantry 2 via the control interface 62. 

20 The data acquisition unit 26, X-ray controller 28, collimator controller 30, and 

rotation controller 36 which are incorporated in the scanner gantry 2 are controlled via the 
control interface 62. Illustrating the connections among these components and the control 
interface 62 will be omitted. 

A data acquisition buffer 64 is also connected to the data processing unit 60. The 

25 data acquisition unit 26 incorporated in the scanner gantry 2 is connected to the data 
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acquisition buffer 64. Data acquired by the data acquisition unit 26 is transferred to the 
data processing unit 60 via the data acquisition buffer 64. 

The data processing unit 60 reconstructs an image using transmitted X-ray signals, 
that is, projection information transferred via the data acquisition buffer 64. A storage 
5 device 66 is also connected to the data processing unit 60. Projection information 
transferred via the data acquisition buffer 64, tomographic image information resulting 
from image reconstruction, and programs that realize the capabilities of the system are 
stored in the storage device 66. 

A display device 68 and an operating device 70 are also connected to the data 

10 processing unit 60. Tomographic image information and other information transferred 
from the data processing unit 60 are displayed on the display device 68. An operator 
handles the operating device 70 so as to enter various instructions or information which is 
then transferred to the data processing unit 60. The operator uses the display device 68 
and operating device 70 to interactively operate the X-ray CT system. Incidentally, the 

15 scanner gantry 2, a radiographic table 4, and the operating console 6 constitute an 
acquisition system that scans a subject or a phantom so as to acquire tomographic image 
data. 

Fig. 2 is a functional block diagram showing a portion of the data processing unit 
60 relating to a beam-hardening post-processing method in accordance with the present 
20 embodiment. 

The data processing unit 60 comprises a beam-hardening correction block 201, a 
first fitting block 202, a second fitting block 204, a correction coefficient modification 
block 205, and an image reconstruction block 206 which deal with projection information 
stored in the storage device 66. 
25 The beam-hardening correction block 201 corrects projection information stored in 
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the storage device 66 in terms of the beam-hardening effect. Assuming that projection 
information values detected on each of the channels of the X-ray detector 24 are Ih and 
data corrected in terms of the beam-hardening effect is IC 5 correcting projection 
information in terms of the beam-hardening effect is expressed as follows: 

5 

IC = B 0 • Ih + Bi • Ih 2 + B 2 • Ih 3 4- B 3 • Ih 4 (1) 

where B 0 to B 3 denote correction coefficients. These correction coefficients are finalized 
in relation to each of the channels of the X-ray detector according to a method described 

10 in, for example, Patent Document 1, and stored in the form of a correction coefficient table 
in the storage device 66. 

The first fitting block 202 smoothes projection information values, which are 
sampled in relation to all the respective views and each channel or which are sampled in 
relation to each view and all the respective channels, from projection information stored in 

15 the storage device 66. The first fitting block 202 averages the projection information 
values that are sampled in relation to each channel and all the respective views or that are 
sampled in relation to all the respective channels and each view. Otherwise, the first 
fitting block 202 fits a high-order function to the projection information values sampled in 
the direction of channels or in the direction of views so as to average the projection 

20 information values. 

A function resulting from the fitting does not reflect high-frequency components of 
acquired signals that are higher than frequency components determined with the order of 
the function. The fitting therefore provides the same effect as smoothing. 

The second fitting block 204 fits a linear or high-order function to projection 

25 information values which are detected on each of the channels of the X-ray detector 24 
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and to which the first function is fitted by the first fitting block 202. This results in the 
same correction coefficients as the ones provided by the expression (1) solved by the 
beam-hardening correction block 201. 

The correction coefficient modification block 205 modifies the correction 
5 coefficient, which the second fitting block 204 has just calculated as a function, using a 
correction coefficient which the second fitting block 204 has already calculated as a 
function using projection information acquired from other phantom, if necessary (if the 
diameter of the phantom is so large that projection information reflects low signal-to-noise 
ratios). A concrete process will be described later. 

10 The image reconstruction block 206 uses a sinogram, which is produced using 

projection information composed of a plurality of views and stored in the storage device 
66, to reconstruct a tomographic image of a subject or a phantom. For image 
reconstruction, for example, the filtered back projection technique is adopted. A 
reconstructed image is displayed on the display device 68. 

15 Next, an example of a procedure for calculating a correction coefficient needed for 

correction that is performed as post-processing in order to cope with the beam-hardening 
effect in the X-ray CT system 1 will be described with reference to the flowchart of Fig. 3. 
Scanning a phantom 

First, a phantom is scanned (step S501). Specifically, a phantom is disposed at a 
20 position off the center of the X-ray field in the bore 29. The phantom is made of a 
polypropylene or the like and shaped like a cylinder. Phantoms having various diameters 
are available. The present embodiment employs phantoms whose diameters are 35 cm and 
48 cm. The phantoms whose diameters are 35 cm and 48 cm are scanned in that order. 

When a phantom is scanned, the phantom is disposed at a position off the center of 
25 the X-ray field in the bore 29. This is because the length of a path traced by X-rays 

10 



16CT03014 

transmitted by the phantom can be differentiated from view to view. In order to highly 
precisely correct projection information values in terms of the beam-hardening effect, the 
largest possible number of different projection information values is needed in relation to 
each channel. When a phantom is disposed at a position off the center of the X-ray field in 
5 the bore 29, numerous different projection information values can be acquired from one 
phantom. 

Fig. 4 shows a phantom 310 disposed in the bore 29 of the scanner gantry 2. 
The phantom 310 has a circular section, and the center of the phantom 310 is 
located at a position different from the center of the X-ray field in the bore 29. An X-ray 
10 fan beam into which X-rays generated by the X-ray tube 20 are recomposed by a bow-tie 
filter is transmitted by the phantom 310 and detected by the X-ray detector 24. 

The X-ray detector 24 has a plurality of X-ray detection elements set in array in a 
direction in which the X-ray fan beam spreads. Projection information acquired from the 
phantom 310 is detected on the channels assigned to the X-ray detection elements set in 
15 array. Herein, the X-ray tube 20 and collimator 22 are opposed to the X-ray detector 24 
with the bore 29 between them. The X-ray tube 20, collimator 22, and X-ray detector 24 
are rotated about the bore 29 without any change in their relative positions while being 
incorporated in the rotary assembly 34, whereby projection information is acquired. 
Projection information is acquired for each view number associated with a rotation angle, 
20 and one sinogram is produced. 

Fig. 5(A) shows an example of a sinogram produced using the phantom 310. The 
sinogram comprises a projection information section defined around the center of one 
dimension of the sinogram associated with channel numbers, and an air data section 
defined on both sides of the projection information section along the dimension associated 
25 with channel numbers. The phantom 310 is disposed off the center of the X-ray field. 
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Therefore, the channel numbers determined with the channel-direction width of the 
projection information section vary with a change in the rotation angle of the rotary 
assembly 34, that is, a change in the view number. As shown in Fig. 5(A), the projection 
information section is tortuous in the direction of view numbers. For the same reason, the 
5 channel-direction width of the projection information section varies with a change in the 
view number. 

Fig. 5(B) is a graph whose axis of abscissas indicates channel numbers and whose 
axis of ordinates indicates projection information values, wherein projection information 
of view number j indicated in Fig. 5(A) is pointed out. The projection information value is 

10 proportional to a length in the phantom 310 over which an X-ray beam is transmitted. X- 
rays passing through near the center of the phantom 3 10 are transmitted over a large length 
and provide large projection information values. X-rays passing through near the 
periphery of the phantom 310 are transmitted over a small length and provide small 
projection information values based on which a semi-circular projection image like the 

15 one shown in Fig. 5(B) is displayed. 

Now, a projection information value indicated with view number j and channel 
number i will be taken for instance. An X-ray beam indicated with a dashed line in Fig. 4 
provides a view of view number j and falls on a channel of channel number i of the X-ray 
detector 24. At this time, a length in the phantom 310 over which the X-ray beam is 

20 transmitted shall be 1. The length 1 and a projection information value h detected on 
channel i and indicated in Fig. 5(B) has the relationship expressed below. 

1 cc h 

25 Referring to Fig. 4, the phantom 310 is located at a position off the center of the X- 
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ray field. The length 1 over which the X-ray beam is transmitted and routed to the channel 
i varies depending on a view. Consequently, the projection information value h detected 
on the channel i and indicated in Fig. 5(B) varies depending on a view. 

Fig. 5(C) is a graph whose axis of abscissas indicates view numbers and whose axis 
5 of ordinates indicates projection information values, wherein a projection information 
value detected on the channel of channel number i indicated in Fig. 5(A) is pointed out. 
The projection information value varies depending on a view number, and is proportional 
to a length in the phantom 310 over which an X-ray beam is transmitted. Therefore, the 
projection information value is provided as a function represented by a curve having 
10 cyclicity as shown in Fig. 5(C). 

Fig. 6 is a block diagram showing files stored in the storage device 66. 

By scanning the phantom as mentioned above, first projection information 601 
shown in Fig. 6 is acquired. 

Preprocessing 

15 Thereafter, a sinogram produced by scanning the phantom is preprocessed (step 

S502). 

The preprocessing is performed on the sinogram produced from the first projection 
information, and includes noise removal and sensitivity correction. 
Correcting data in terms of the beam-hardening effect 
20 Thereafter, using the expression (1), the projection information values Ih are 

corrected in terms of the beam-hardening effect in order to produce corrected projection 
information values Ic (step S503). 

Consequently, second projection information 602 shown in Fig. 6 is produced. 
This file has the beam-hardening effect generally removed from it. However, the beam- 
25 hardening effect slightly remains because of a difference of each channel of the X-ray 
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detector 24 from the others. 

Fig. 7(A) illustratively shows an example of the second projection information 602. 
The second projection information 602 is plotted generally semi-circularly because it is 
acquired from a phantom having a circular section. The projection information values Ic 
5 sampled in relation to some channels are plotted as pulse-like fluctuations because of a 
difference in the sensitivity to X-rays of each channel from the others. This is a 
phenomenon inherent to a channel, and must therefore be corrected in relation to each 
channel. 

Fig. 8(A) illustratively shows an example of projection information values sampled 
10 from the second projection information 602 in the direction of views in relation to one 
channel. The projection information values Ic sampled in relation to some views are 
plotted as pulse-like fluctuations. 

Smoothing in the direction of channels 

Thereafter, the first fitting block 202 smoothes the second projection information 
15 602 in the direction of channels (step S504). 

Consequently, third* projection information 603 shown in Fig. 6 is produced. The 
third' projection information 603 has the projection information values Ic, which are 
plotted as fluctuations deriving from a difference of one channel from the others, 
smoothed and removed. Fig. 7(B) illustratively shows an example of the third' projection 
20 information 603. A semi-circular curve alone is plotted with the projection information 
acquired from a phantom having a circular section. 
Smoothing in the direction of views 

Thereafter, the first fitting block 202 smoothes the third' projection information 603 
in the direction of views (step S505). 
25 Consequently, third projection information 604 shown in Fig. 6 is produced. The 
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third projection information 604 has projection information values, which are plotted as 
fluctuations deriving from a difference of one view from the others detected on one 
channel, smoothed. 

Fig. 8(B) illustratively shows an example of the third projection information. The 
5 projection information values that are sampled in relation to one channel in the direction 
of views and that exhibit cyclicity are smoothed. 
Calculating a correction coefficient 

Thereafter, the second fitting means 204 calculates a correction coefficient using 
the second projection information 602 and third projection information 603 (step S506). 

10 Herein, projection information values sampled from the second projection 

information in relation to channel number i shall be S(j) and projection information values 
sampled from the third projection information in relation thereto shall be F(j). Fig. 9(A) 
indicates the projection information values sampled in relation to all the view numbers, 
wherein the axis of abscissas indicates the projection information values S(j) and the axis 

15 of ordinates indicates the projection information values F(j). The projection information 
values are generally plotted alone a straight line that passes an origin. The straight line 
represents a correction function relevant to the channel i. Assuming that the slope of the 
straight line is Kji, Kji is regarded as correction coefficient information 605 shown in Fig. 
6. Incidentally, S(j) and F(j) have the relationship expressed below. 

20 F(j)/S(j) = Ki 

The projection information values Ic sampled in relation to channel i from 
projection information that is acquired from a subject and that is corrected in terms of the 
beam-hardening effect are multiplied by the correction coefficient Ki of a linear function. 
Consequently, projection information values Ip are produced as projection information 

25 values that have been smoothed and corrected after being acquired from the subject. 
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The projection information values S(j) are proportional to the length 1 over which 
X-rays are transmitted and which is indicated in Fig. 4. Therefore, the projection 
information values S(j) and the adopted domain depend on the diameter of the phantom 
310 and the position of the phantom in the bore 29. 
5 Thereafter, the correction coefficient Ki is preserved as correction coefficient 

information 605 in the storage device 66 (step S507). 

Thereafter, it is determined whether a phantom having a different size is scanned in 
order to acquire data (step S508). 

For example, after a phantom having a diameter of 35 cm is scanned in order to 
10 produce correction coefficient information, if a phantom having a diameter of 48 cm 
should be scanned in order to produce correction coefficient information, the phantom 
having a diameter of 48 cm is disposed in the X-ray CT system 1 . The foregoing steps 
S501 to S507 are then resumed. 

Fig. 10 is a graph indicating an example of the correction coefficient information 
15 605 produced in relation to each channel. Fig. 10(a) indicates correction coefficients 
calculated using the phantom having a diameter of 35 cm, and Fig. 10(b) indicates 
correction coefficients calculated using the phantom having a diameter of 48 cm. 

The correction coefficients calculated using the phantoms whose diameters are 35 
cm and 48 cm are plotted as similar characteristic curves in relation to the channels. 
20 However, the correction coefficients calculated using the phantom having a diameter of 48 
cm assume larger values. 

The X-ray detector 24 comprises a plurality of detection modules each assigned a 
predetermined number of channels. The detection modules are manufactured as united 
bodies in units of sixteen channels. In other words, a predetermined number of detection 
25 modules each assigned sixteen channels is arrayed in order to construct the X-ray detector 
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24. 

Therefore, the values detected on the channels within the same detection module 
reflect similar detection characteristics. In other words, the detection modules exhibit 
different detection characteristics. 
5 Due to the above characteristics of the X-ray detector 24, the correction coefficient 

values arranged orderly in relation to the channels as shown in Fig. 10(a) and Fig. 10(b) 
contain the dependencies on the detection characteristics of the modules. 

Moreover, as seen from Fig. 10(a) and Fig. 10(b), the correction coefficient values 
arranged orderly in relation to the channels contain the reflections of high-frequency 
10 components that depend on the detection characteristics of the respective channels. 

The larger the diameter of a phantom, the lower the signal-to-noise ratios reflected 
in projection information detected by the X-ray detector 24. Besides, the precision in 
projection information detected on each channel gets lower. The reflections of high- 
frequency components in the correction coefficient values indicated in Fig. 10(b) include 
15 more reflections of noises than the reflections thereof in the correction coefficient values 
indicated in Fig. 10(a). 

Therefore, when the correction coefficients calculated using, for example, the 
phantom having a diameter of 48 cm are adopted as they are, a reconstructed image may 
suffer a defect such as a ring artifact. 
20 Modifying a correction coefficient 

After the correction coefficient information 605 is produced using the phantom 
having a diameter of 48 cm by performing the steps S501 to S507, it is judged whether the 
correction coefficient information 605 should be modified (step S509). 

At this time, the correction coefficient information 605 reflecting low signal-to- 
25 noise ratios as indicated in Fig. 10(b) is modified. Correction coefficient information 
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produced using the phantom having a small diameter of 35 cm is not modified. 

Fig. 11 includes graphs for explaining an example of a method of modifying 
correction coefficient information 605. 

Referring to Fig. 11(a), SA indicates correction coefficients obtained by removing 
5 the reflections of high-frequency components from the correction coefficients that are 
calculated using the phantom having a diameter of 48 cm and indicated in Fig. 10(b). 

The correction coefficients SA contain the dependencies on the detection 
characteristics of the detection modules but do not contain the reflections of high- 
frequency components that depend on the detection characteristics of the channels and that 
10 degrade signal-to-noise ratios. 

Referring to Fig. 11(b), a wave SB indicates the reflections of high-frequency 
signal components separated from the correction coefficients that are calculated using the 
phantom having a diameter of 35 cm and indicated in Fig. 10(a). 

The wave SB indicates the reflections of high-frequency signal components in the 
15 correction coefficients that are calculated using the phantom having a diameter of 35 cm 
and that depend on the detection characteristics of the respective channels. The reflections 
of high-frequency components contain the reflections of noises at a smaller ratio than the 
reflections of high-frequency components in the correction coefficients calculated using 
the phantom having a diameter of 48 cm. 
20 During modification of correction coefficient information 605 employed in the 

present embodiment, as indicated in Fig. 11(c), the correction coefficients SA from which 
the reflections of high-frequency components are removed are synthesized with the 
reflections of high-frequency components in the correction coefficients calculated using 
the phantom having a diameter of 35 cm. The resultant correction coefficients are adopted 
25 as new correction coefficients SC. 
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This results in the improved precision in correction coefficient information 
calculated using the phantom having a diameter of 48 cm. As shown in Fig. 6, the 
modified correction coefficient information is stored as modified correction coefficient 
information 606 in the storage device 66. 
5 Fig. 12 shows an example of correction coefficients calculated using two phantoms 

having different diameters (35 cm and 48 cm). 

Referring to Fig. 12, domain A comprises projection information values acquired 
from the phantom having a diameter of 35 cm, and domain B comprises projection 
information values acquired from the phantom having a diameter of 48 cm. 
10 Fitting a high-order function 

Thereafter, a high-order function is fitted to a correction coefficient calculated 
using the phantom having a diameter of 35 cm and a correction coefficient calculated 
using the phantom having a diameter of 48 cm (step S51 1). 

Fig. 13 shows an example of fitting of a high-order function performed using the 
15 correction coefficient information indicated in Fig. 12. 

A third-order function given as expression (2) below is fitted to the correction 
coefficient value calculated from domain A and the correction coefficient value calculated 
from domain B, whereby correction coefficients K0, Kl, and K2 are determined. 

20 If=K0- SG) + K1- S(j) 2 + K2- SG) 3 (2) 

Incidentally, the correction coefficient calculated from domain A composed of 
small projection information values is thought to be more highly precise than the 
correction coefficient calculated from domain B composed of large projection information 
25 values. Therefore, weights may be applied to the respective domains so that the third- 
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order function will be more highly precisely fitted to the correction coefficient calculated 
from the domain A, and the correction coefficients in the expression (2) may then be 
determined. 

Thereafter, as shown in Fig. 6, high-order correction function information 607 
containing the correction coefficients K0, Kl, and K2 is preserved in the storage device 66 
(step S509). The process is then terminated. 

When a subject is scanned, projection information values Ic acquired from the 
subject and corrected in terms of the beam-hardening effect are corrected using correction 
coefficients K0, Kl, and K2 that are calculated in relation to each channel. Corrected 
projection information values If are then calculated according to the expression (2). Based 
on the projection information values If, the image reconstruction block 206 reconstructs an 
image so as to acquire tomographic image information. 

As mentioned above, according to the present embodiment, correction coefficient 
information that is calculated using a large phantom and that reflects low signal-to-noise 
ratios is corrected using correction coefficient information that is calculated using a 
phantom smaller than the phantom and that reflects high signal-to-noise ratios. 
Consequently, projection information values corrected in terms of the beam-hardening 
effect and sampled in relation to each channel can be corrected highly precisely. 

Many widely different embodiments of the invention may be configured without 
departing from the spirit and the scope of the present invention. It should be understood 
that the present invention is not limited to the specific embodiments described in the 
specification, except as defined in the appended claims. 
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